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Hamiltonian perturbations of the simplest hyperbolic equation ut + a{u)u x = 
are studied. We argue that the behaviour of solutions to the perturbed equation 
in . near the point of gradient catastrophe of the unperturbed one should be essen- 

tially independent on the choice of generic perturbation neither on the choice of 



generic solution. Moreover, this behaviour is described by a special solution to 
an integrable fourth order ODE. 



1 Introduction 



In the present work we continue the study of Hamiltonian perturbations of hyperbolic 
PDEs initiated by the paper [TU]. We consider here the simplest case of a single 
equation in one spatial dimension 

Ut + a(u)u x + e [bi(u)u xx + b 2 {u)u x ] + e 2 [b 3 (u)u xxx + b A (u)u x u xx + b 5 (u)u x ] H = 0. 

(1.1) 

Here e is a small parameter; the coefficient of e k is a graded homogeneous polynomial 
in the derivatives u x , u xx , ... of the total degree (k + 1), 

deg = 72,, n > 0. 

The unperturbed equation 

u t + a{u)u x = (1.2) 

can be considered as the simplest example of a nonlinear hyperbolic system; the smooth 
functions bi(u), b 2 (u) etc. determine the structure of the perturbation. 
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Such expansions arise, e.g., in the study of the long wave (also called dispersionless) 
approximations of evolutionary PDEs; see section El below for other mechanisms that 
yield perturbed equations of the form (jl.lj) . 

The unperturbed equation (jl.2j) admits a Hamiltonian description of the form 

5H 

u t + {u(x),H } = u t + d x ° = (1.3) 

ou(x) 



H = J f(u) dx, f"(u) = a(u) 

{u{x),u{y)} = 5\x-y) (1.4) 

The perturbed equations of the form (jl.lj) are considered up to equivalencies defined 
by Miura-type transformations |H] of the form 

u^u + J2e k F k (u;u x ,...,u (k} ) (1.5) 

k>l 

where F k (u;u x , . . . ,u^) is a graded homogeneous polynomial in the derivatives u x , 
u xx , ... of the degree 

deg F k = k. 

Using results of [T3] (see also [HHH]) one can show that any Hamiltonian perturbation 
of the equation (jl.2jl can be reduced to the form 

v>t + d x J^- = 0, H = H + eH 1 + e 2 H 2 + ... 

(1.6) 

H k = J h k (u; u x ,..., ii (fc) ) dx, deg h k (u; u x , . . . , u {k) ) = k. 

Recall that for H = J h(u; u x , u xx , . . . ) dx 

SH 

Su(x) 

where 



Eh 



du du x x du xx 

is the Euler - Lagrange operator. The following well known property of the Eu- 
ler - Lagrange operator will be often used in this paper: £ h = iff there exists 
hi = hi(u; u x , . . . ) such that h = const + d x h\. Note that we do not specify here the 
class of functions u(x). The Hamiltonians H = H[u] can be ill defined (e.g., a diver- 
gent integral) but the evolutionary PDE (jl.6j) makes sense. The crucial point for the 
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subsequent considerations is the following statement (see, e.g., [Zj): for two commuting 
Hamiltonians 

SH „ 5F 



the evolutionary PDEs 

w * + ^77~\ = and u s + Ox^-rr = 
ow(a;j on(x) 

commute, 

(U t )s = (U s ) t . 

For sufficiently small e one expects to see no major differences in the behaviour 
of solutions to the perturbed and unperturbed equations (jl.lj) and (jl.2j) within the 
regions where the ^-derivatives are bounded. However the differences become quite 
serious near the critical point (also called the point of gradient catastrophe) where the 
derivatives of solution to the unperturbed equation tend to infinity. 

Although the case of small viscosity perturbations has been well studied and un- 
derstood (see [3] and references therein), the critical behaviour of solutions to general 
conservative perturbations to our best knowledge has not been investigated (see 
the papers [T21 EE7J CHI EH 123 122 I2H1 123 for the study of various particular cases). 

The main goal of this paper is to formulate the Universality Conjecture about the 
behaviour of a generic solution to the general perturbed Hamiltonian equation near the 
point of gradient catastrophe of the unperturbed solution. We argue that, up to shifts, 
Galilean transformations and rescalings this behaviour essentially does not depend on 
the choice of solution neither on the choice of the equation (provided certain genericity 
assumptions hold valid). Moreover, this behaviour near the point (xQ,t ,uo) is given 

by 



u 



e 2 / 7 U (fee" 6 / 7 (x — a (t — t ) — x ); ce 4//7 (£ — t )) + O (e 4 / 7 ) (1.7) 



where U = U(X; T) is the unique real smooth for all X G K solution to the fourth 
order ODE 



X = TU- 



-U 3 + — (U' 2 + 2U U") + — U IV 
6 24 v ; 240 



V=% etc. (1.8) 



depending on the parameter T. Here a, b, c are some constants that depend on the 
choice of the equation and the solution, a = a(vo). 

The equation ()1.8|) appeared in |lj (for the particular value of the parameter T = 0) 
in the study of the double scaling limit for the matrix model with the multicritical index 
m = 3. It was observed that generic solutions to (jl.8|) blow up at some point of real line; 
the conjecture about existence of a unique smooth solution has been formulated. To 
our best knowledge, this conjecture remains open, although there are some supporting 
evidences 
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The present paper is organized as follows. In Section 2 we classify all Hamiltonian 
perturbations up to the order e 4 . They are parametrized by two arbitrary functions 
c(u), p(u). For the simplest example the perturbations of the Riemann wave equation 
u t + u u x = read 

(1.9) 

+2p'(5u xx u xxx + 3u x u xxxx ) + p"{7u x u 2 xx + <ou 2 x u xxx ) + 2p"'u 3 x u xx ] = 0. 

For c(u) = const, p(u) = this is nothing but the Korteweg - de Vries (KdV) equation; 
for other choices of the functions c(u), p(u) it seems not to be an integrable PDE. 
Remarkably, for arbitrary choice of the functional parameters the perturbed equation 
possesses an infinite family of approximate symmetries (see [2J HH 1221 EH] for discussion 
of approximate symmetries). In principle our approach can be applied to classifying 
the Hamiltonian perturbations of higher orders. However, higher order terms do not 
affect the type of critical behaviour. 

In Section 3 we establish an important property of quasitriviality of all perturbations 
(cf. |Hl EH I2Z1 ) • The quasitriviality is given by a substitution 

u h-» u + e 2 K 2 {u\ u x , u xx , u xxx ) + e 4 K 4: (u; u x , ... , w (6) ) (1.10) 

that transforms, modulo 0(e 6 ) the unperturbed equation (J1.2)) to (jl.6|) . Here the 
functions K 2 and K4 depend rationally on the ^-derivatives. We also formulate the 
first part of our Main Conjecture that says that, for sufficiently small e the solution 
to the perturbed system exists at least on the same domain of the (x, t)-plane where 
the unperturbed solution is defined. In Section 4 we briefly discuss existence of a 
bihamiltonian structure compatible with the perturbation (see also Appendix below). 
Some examples of perturbed Hamiltonian equations are described in Section 5. In 
Section 6 we recollect some properties of the ODE ()1.8|) . Finally, in Section 7 we 
give the precise formulation of the Universality Conjecture and give some evidences 
supporting it. Because of lack of space we do not consider the numerical evidences 
supporting the idea of Universality; they will be given in a subsequent publication (see 
also ^H])- In the last section we outline the programme of further researches towards 
understanding of universality phenomena of critical behaviour in general Hamiltonian 
perturbations of hyperbolic systems. 

Acknowledgments. This work is partially supported by European Science Foun- 
dation Programme "Methods of Integrable Systems, Geometry, Applied Mathematics" 
(MISGAM), Marie Curie RTN "European Network in Geometry, Mathematical Physics 
and Applications" (ENIGMA), and by Italian Ministry of Universities and Researches 
(MIUR) research grant PRIN 2004 "Geometric methods in the theory of nonlinear 
waves and their applications" . 
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2 Hamiltonian perturbations of the Riemann wave 
equation 

Let us start with the simplest case of Hamiltonian perturbations of the equation 

v t + vv x = & v t + {v(x),H } = (2.1) 

{v(x),v(y)} = S'(x - y) 

fv 3 , 
H = / — dx 



Lemma 2.1 Up to the order 0(e A ), all Hamiltonian perturbations of \2.1\) can be re- 
duced to the form 

„ 5H 

Ut + Ox——^- = 

du(x) 



H 



cix (2.2) 



where c(u), p(u), s(u) are arbitrary functions. Moreover, the function s(u) can be 
eliminated by a Miura-type transform. 

Proof The Hamiltonian must have the form 

H = H + eH 1 + --- + e 4 # 4 

where the density of H k is a graded homogeneous polynomial of the degree k. So, the 
density of Hi is a total derivative: 

Hi = a(u)u x dx, a(u)u x = d x A{u), A'(u) = a(u). 



The density of the Hamiltonian H 2 modulo total derivatives must have the form 



c(u) 



24 x 

for some function c(u). Similarly, H3 must have the form 

H > = J 

Here C\{u) is another arbitrary function. 



Let us show that H3 can be eliminated by a Miura-type transform. Let us look for 
it in the form 

u 1 ► u + e{u(x),F} + — {{u(x),F}, F} + . . . (2.3) 

choosing 

F = e 2 J a{u)u 2 x dx. 

Such a transformation preserves the Poisson bracket. The change of the Hamiltonian 
H will be given by 

5H = e{F,H} + 0(e 4 ). 

At the order e 3 one has 



6H = e 3 



-a'(u)ul - d x (au x ) 



e 3 



uu x dx = — I a(u)u x dx. 



So, choosing a(u) = —2c\{u) we kill the terms cubic in e. 

The rest of the proof is obvious: in order 4 all the Hamiltonians have the form 



for some functions p(u), s(u). The last term can be killed by the canonical transfor- 
mation of the form (J2.3)) generated by the Hamiltonian 



e 3 

-— / s(u)uldx. 



The lemma is proved. 

Choosing s(u) = one obtains the family (|1.9|) of Hamiltonian perturbations of the 
Riemann wave equation depending on two arbitrary functions c = c(u) , p = p(u) . 

We will now compare the symmetries of (j2.1J) and those of the perturbed system 
()2.2|) . It is easy to see that the Hamiltonian equation 

v s + a(v)v x = <£> v s + {v(x),H° f } = (2.4) 
H° f = J f(v)dx, f"(v) = a(v) 

is a symmetry of (|2.1j) for any a(v), 

{v t ) s = (v s )t- 

Moreover, the Hamiltonians H® commute pairwise, 

{H° f ,H° g } = V/ = /(«), Vg = g(u). 
This family of commuting Hamiltonians is complete in the following sense. 
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Lemma 2.2 The family of commuting Hamiltonians H® is maximal, i.e., if H = 
J h(u; u x , u xx , . . . ) dx commutes with all functionals of the form H® then 

h(u; u x , u xx , ■■■) = g{u) + d x {. . . ) 

for some function g(u) . 

We will now construct a perturbation of the Hamiltonians Hf preserving the com- 
mutativity modulo 0(e 6 ). Like in Lemma 12. II one can easily check that all the pertur- 
bations up to the order 4 must have the form 

Hf = J {/(«) - e^J^ul + e 4 [p f (u)ul x + *,(«)<] } dx 

for some functions Cf(u), Pf(u), Sf(u). To ensure commutativity one has to choose 
these functions as follows. 

Lemma 2.3 For any f = f{u) the Hamiltonian flow 



5H 



"Hi 2 , 4 



/(4) 



cc"/ (4) + cc'fV + If® + pJ^ + pfV _ sfn 



1152 1152 3456 6 6 

is a symmetry, modulo 0(e 6 ), of Moreover, the Hamiltonians Hf commute 

pairwise: 

{H f ,H 9 } = 0(e ti ) 
for arbitrary two functions f{u) and g{u). 

Proof One has to check the identity 

S (^d x ^A = 

where £ is the Euler - Lagrange operator. We leave this calculation as an exercise for 
the reader. 

Observe that for / = ^ the Hamiltonian Hf coincides with (J2.2|) . Also for / = u 
(the Casimir of the Poisson bracket) and f = \ (the momentum) the perturbation is 
trivial, 

Hf = H° f . 

We do not know under what conditions on the functional parameters c(u), p(u) 
higher order perturbations can be added to the Hamiltonians (J2.5)) preserving the 
commutativity. The examples of Section show that this can be done at least for 
some particular choices of the functions. However, the remark at the end of Section 0] 
suggests that the answer is not always affirmative. 
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3 Solutions to the perturbed equations. Quasitriv- 
iality 



Next question: existence of solutions to the perturbed equation for t < tc- We will 
construct a formal asymptotic solution to (|2.2|) (and also to all commuting flows (|2.5j) ) 
valid on the entire interval t < tc- The basic idea: find a substitution 

v u = v + 0(e) 

that transforms all solutions to all unperturbed equations of the form (|2.4J) to solutions 
to the corresponding perturbed equations (|2.5j) . 

Quasitriviality Theorem There exists a transformation 



u 



Y,£ k F k (u; u x ,...,u^), 



(3.1) 



fc=i 



where Fk are rational functions in the derivatives homogeneous of the degree k, in- 
dependent of f = f{u), that transforms all monotone solutions of \2-4\j to solutions, 
modulo 0(e 6 ) ; of Q2.ty) and vice versa. 

The general quasitriviality theorem for evolutionary PDEs admitting a bihamilto- 
nian description was obtained in ^Uj 1 . As we do not assume a priori existence of a 
bihamiltonian structure (see, however, the next section), we will give here a direct 
proof of quasitriviality for the family of commuting Hamiltonians (|2.5)1 . 

For convenience we chose 

c(u) c"'(u) 



s(u) 



3456 



Theorem 3.1 Introduce the following Hamiltonian 
K -- 

Then the canonical transformation 



—ec(u) u x \ogu x + -J^-L^EL _ ^__L^_ 
24 K J x B I 5760 ul 4 u x 



dx. 



satisfies 



u i— > v = u + e{u(x), K} H {{u(x), K}, K} 



H f = J f(v)dx + 0(e 6 ) V/(«). 



-"^In a very recent paper |27j the quasitriviality result was proved, in all orders in e, for an arbitrary 
perturbation of the Riemann wave equation Vt + v v x — 0. It has also been shown that the same 
transformation trivializes also all symmetries of the perturbed equation. 
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The inverse transformation is the needed quasitriviality. It is generated by the 
Hamiltonian 



-K 



24 V 5760 vt 4 v 



that is 



v i — > u = v — e{v(x),K} + — {{v(x), K}, K} + 

= v + —,d x ( c-^- + c'v x | + e^d x 
24 V v r 



7 U ^ X X X 



+ 



360 v$ 1920 vl 1152 v*, a 



+CC 



47 f. 



37 

^ XX ^ XXX ^ XXXX \ /2 ^ ^XX'i 

+ — — — + c 



5760 w, 3 2880 7V 2 1152 



384 5760 v,, 



+ cc" 



v. 



144 360 u, 



+ iT52 ( 7 °' C " Vx Vxx + V * + 6 ° C '" ^ ^ + C ' C "' Vx> + C ° (4) 1,2:3 
/ ?;„, 3 



^ XX ^ XXX ^ xxxx 

+ — I + P'v xxx + p" 



2v x 3 v x 2 2v x 

In this formula c = c(v), p = p(v). 

Main Conjecture, Part 1. Let v = v(x, t) be a smooth solution to the unperturbed 
equation Vt + a(v) v x — defined for all x £l and < t < to monotone in x for any 
t. Then there exists a solution u = u(x,t; e) to the perturbed equation 

SHj 



u t + d x 



Su(x) 



0, f"{u) = a(u) 



defined on the same domain in the (x,t) -plane with the asymptotic at e — ► of the 
form (TOj) . 



4 Are all Hamiltonian perturbations also bihamil- 
tonian? 

All unperturbed equations 

v s + a(v) v x = 

are bihamiltonian w.r.t. the Poisson pencil (see the definition in jUJ) 

{v{x),v(y)} 1 = 5'{x - y), {v{x),v{y)} 2 = q(v(x))S'(x - y) + -q'{v)vj{x - y) (4.1) 
for an arbitrary function q(u), 

v s + {v(x),H 1 } 1 = v s + {v(x),H 2 } 2 = 0, Hi = J fx(v) dx, H 2 = J f 2 {v) dx 

f';{v) = a{v) = q{v)f' 2 \v)+ l -q\v)f 2 {v). 



(3.2) 
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To show that (|4.1J) is a Poisson pencil it suffices to observe that the linear combination 



{v(x),v(y)} 2 - A {v(x),v(y)} 1 = (q(v{x)) - A) S'(x - y) + -q'{v)v x 5(x - y) (4.2) 



is the Poisson bracket associated jTTj with the flat metric 



ds 1 



dv 2 



q(v) - A' 



Theorem 4.1 For c(u) ^ the commuting Hamiltonians \2. ,5)) admit a unique bi- 
hamiltonian structure obtained by a deformation of fl^.i| ) with q(u) satisfying 



p(u) 



960 



c g 
5 

c g' 



s hi 



(4.3) 



The proof of this result along with the explicit formula for the deformed bihamil- 
tonian structure is sketched in the Appendix below. 

The assumption c ^ is essential: one can check that for c(u) = the Hamiltonians 
(J2.5J) commute, modulo 0(e 6 ), only w.r.t. the standard Poisson bracket (jl.4j) . On the 
other side it turns out that for this particular choice of the functional parameters the 
deformation of commuting Hamiltonians cannot be extended to the order 0(e 8 ). 



5 Examples 



Example 1. For c(u) = cq = const, p(u) = s(u) = one obtains from ()2.2|) the KdV 
equation 



u t + u u x + c —u xxx = 0. 



Choosing in (|2.5|) 



/(«) 



fc+2 



(fc + 2)! 

one obtains the Hamiltonians of the KdV hierarchy 



— + d — 

x Su(x) 



0, A* 



/ifc k > 



it 



fc+2 



- c 



e 2 -u fc 1 



u 



fc-2 



5 (fc-2) 



(fc + 2)! u 24(fc-l)!"* ' ~ u 96 
The quasitriviality transformation ()3.2|) takes the form [2 

7 " :i 7v xx v 



U xx 



a 



fc-4 



36(fc-4)! 



+ 0(e 6 



w i— ► u — v + 9„ 



v 3 

u xx 



— c log v x + c 2 e 4 . 

24 6 \360v$ 1920 u2 1152 w 2 



+ 0(e 6 ). (5.1) 
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Example 2. The Volterra lattice 

q n = qn(q n +i - g n -i) (5-2) 

(also called difference KdV) has the following bihamiltonian structure [T3] 

) (5-3) 



q n = {q n ,H 1 } 1 , H 1 = ^^2\ogq n 



{q n , q m }2 



(<5n,m+l — <^n,m-l) + 2^»Vn+2 — T^™- 7 ™ -2 H^) 



g n = {q n , H 2 } 2 , H 2 = y^ j q n 
After substitution 

q n = e" 

and division by 4e one arrives at the following bihamiltonian structure 



1 / e 2 

= ^ [<J(x - y + e) - 8(x - y - e)] = 5 (x - y) + —8 (x - y) + . . . (5.5) 



{v(x),v(y)} 2 = (1 - e"W) 8\x - y) - \e v v x 5{x - y) 



(5.6) 



L(2-5e v )8'"(x-y)-^vJ"(x-y) 



3 1 

-~e v (v xx + vl)8'{x -y)- —e v (v xxx + 3v x v xx + vl)8(x - y) 



+ 0(6 4 ). 



To compare this bihamiltonian structure with the one obtained in Theorem 14.11 the 
Poisson bracket ()5.5j) must be reduced to the standard form 



{u(x),u(y} 1 = 8'{x - y) 



(5.7) 



by means of the transformation 



u 



— - — — v = v v xx H v xxxx + U(e ). 

sinhe&. 12 160 V ; 



After the transformation the second bracket takes the form 



{u(x),u(y)} 2 = (1 - e u ^) S'(x - y) - -e u u x 5(x - y) 



(51 



_ € 2 e u(x) 



1 



1.3 1 

1 5 "'{ x -y) + gUx8"(x -y) + y^ Uxx + 5u l) s> ( x ~ v) 



+^(2u xxx + Au x u xx + ul)5(x - y) 



0(e 4 
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We leave as an exercise for the reader to compute the terms of order e 4 and to verify 
that the Poisson bracket (|5.8jl is associated with the functional parameters chosen as 
follows | | 

c(u)=2, P( u ) = ~^ q(u) = l-e u , s(u) = 
Example 3. The Camassa - Holm equation jS] (see also [Ti] ) 



v t ~ e v xx t = -v v x — e 2 
2 



1 

2 



(5.9) 



admits a bihamiltonian description (cf. [21]) after doing the following Miura-type 
transformation 

u = v — e 2 v xx . (5.10) 

The bihamiltonian structure reads 

{u(x),u(y)} 1 = S'(x -y)- e 2 8"'{x - y) (5.11) 



{u{x),u{y)} 2 = u(x)5'(x -y) + -uj{x - y). 



(5.12) 



The Casimir H-\ of the first Poisson bracket analytic in e has the form 



Applying the bihamiltonian recursion procedure one obtains a sequence of commuting 
Hamiltonians = J h^dx of the hierarchy, 

^o = J", h x = hv 3 + uv 2 }, . . . 

I o 



The corresponding Hamiltonian flows 



u tk = {u(x),H k } l = (l-e 2 d 2 )d x 



8u(x) 



read 



u to = u x , u tl = -vv x - e 



2 



The last equation reduces to (|5.9|) after the substitution (|5.10|) . 

To compare the commuting Hamiltonians with those given in ()2.5|) one must first 
reduce the first Poisson bracket to the standard form {u(x), u(y)}i = 8'{x — y) by the 
transformation 



u = (l - e 2 d 2 x ) 1/2 u = u + ^e 2 u xx + ^e 4 u xxxx + 
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After the transformation the Camassa - Holm equation will read 

~^uu x -\- £ (2u x U xx -\- UU XXX ) -\- 6 (5 u xx u xxx ~\~ 3 u x u xxxx ~\~ u u xxxxx ) ~\~ ■ ■ ■ ■ 

It is easy to see that the commuting Hamiltonians of Camassa - Holm hierarchy are 
obtained from ()2.5|) by the specialization 

u 

c(u) =8u, p(u) = — , q(u) = u, s(u) = 0. 
o 

6 Introducing a special function 

Let us remind some properties of the differential equation 
X = TU 



-U 3 + — {U' 2 + 2U U") + — U IV 
6 24 v ; 240 



(6.1) 



often considered as a 4th order analogue of the classical Painleve-I equation. First, 
it can be interpreted as a monodromy preserving deformation of the following linear 
differential operator with polynomial coefficients 

^ = W^ (6.2) 

where the matrix W reads 

, / 12UU' + 8zU' + U'" 2(lQz 2 + 8zU + QU 2 + U" - 60T) 

W 



120 

where 



2% -12 UU' - 8zU' - U- 



ni 



w 21 = 32z 3 - 16z 2 U - 2z{2U 2 + U" + 60 T) + 8U 3 + 2U"U - U' 2 + 120X 
Indeed, it coincides with the compatibility conditions 

W x -U z +[W,U] = 

of the linear system (j6.2j) with 

!t = UV, U=( " ' j (0.3) 



dX 



2U-2z 



Moreover, the dependence of ()6.2j) on T is isomonodromic iff 'the function U(X) depends 
also on the parameter T according to the KdV equation 

U T + UU' + ^-U"' = 0. (6.4) 
13 



This is the spelling of the compatibility condition of the linear system ()6.2|) . ()6.3)1 with 



dT V ' 6 



U' 2U + Az 

lz 2 - AzU - AU 2 - U" -U' 



(6.5) 



The Painleve property readily follows from the isomonodromicity: singularities in the 
complex (X, T)-plane of general solution to (j6.1J) . ()6.4|) are poles |2T)| . 

Main Conjecture, Part 2. T/ie ODE \6.1\) has unique solution U = U (X; T) 
smooth for all real XgI /or a// real values of the parameter T . 

Note that, due to the uniqueness the solution in question satisfies the KdV equation 
fl6~4~j) . 

For T << the solution of interest is very close to the unique root of the cubic 
equation 

U 3 



X~TU , 

6 ' 



that is, 



U 



-T)- 7 w 



w + { _ T y7/2_ 



3w 2 



3 (w 2 + 2) 4 
189u; 4 - 972w 2 + 436 



X = -(-Tf 2 



9 (w 2 + 2) 9 

1 . 
w H — w 
6 



+ 0((-T)- 21 / 2 ) 



(6.6) 



Same is true for any T for \X\ >> 0. For T >> the solution develops oscillations 
typical for dispersive waves within a region around the origin; one can use Whitham 
method to approximate U(X; T) by modulated elliptic functions within the oscillatory 
zone jIHl EH|- Thus the solution in question interpolates between the two types of 
asymptotic behaviour (cf. j2Bj where the role of the special solution U(X; T) in the 
KdV theory was discussed). 

The solutions to the fourth order ODE (J6.1)) can be parametrized [20 by the mon- 
odromy data (i.e., the collection of Stokes multipliers) of the linear differential op- 
erator (j6.3|) with coefficients polynomial in z. The solution corresponding to given 
Stokes multipliers can be reconstructed by solving certain Riemann - Hilbert problem. 
The particular values of the Stokes multipliers associated with the smooth solution in 
question have been conjectured in [20j. 
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7 Local Galilean symmetry and critical behaviour 



We will now proceed to discussing the universality problem. Consider the perturbed 
PDE 

u t + {u(x),H f } = u t + a(u)u x + 0(e 2 ) = 0, f'\u)=a{u). (7.1) 

Let us apply the transformation (jH.2j) to the unperturbed solution v = v(x,t) of 

v t + a(v)v x = (7.2) 

obtained by the method of characteristics: 

x = a(v)t + b(v) (7.3) 

for some smooth function b(v). Let the solution arrive at the point of gradient catas- 
trophe for some x = xq, t = to, v = vq. At this point one has 

x = a(v )t + b(v ) 

= a'(v )t + b'(v ) (7.4) 
= a"(v )t + b"(v ) 

(inflection point). Let us assume the following genericity assumption 

K:=-(a'"(v )t + b"'(v ))^0. (7.5) 

Let us first remind the universality property for the critical behaviour of the un- 
perturbed solutions: up to shifts, Galilean transformations and rescalings a generic 
solution to (|7.2|) near (xq, to) behaves like the cubic root function. We will present this 
well known statement in the following form. Introduce the new variables 

x = x — a (t — t ) — x 
i=t-t 

V = V — Vo. 

Let us do the following scaling transformation 

t^\H (7.6) 



v i— > A 3 v 

Lemma 7.1 After the reseating \7.dj) any generic solution to \7.2j) at the limit A — > 
for t < to goes to the solution of the cubic equation 

x = a' vt — k— . (7.7) 
6 
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In these formulae ao = a(i>o), a' = a'(vo). Note that the inequality 

Ka' >0 (7.8) 

must hold true in order to have the solution well defined for t < t near the point of 
generic gradient catastrophe (|7.4jl . 

To prove the lemma it suffices to observe that, after the rescaling (|7.fij) and division 
by A the equation (|7.H|) yields 

_ g 

x = a' vi- k — + O (A 1/3 ) . 



The parameter k can be eliminated from ()7.7j) by a rescaling. The resulting cubic 
function can be interpreted as the universal unfolding of the A2 singularity p. Our 
basic observation we are going to explain now is that, after a Hamiltonian perturbation 
the A 2 singularity transforms to the special solution of ()1.8|) described above. 

Let us look for a solution to the perturbed PDE (j7.1|) in the form of a formal power 
series 

u = u(x, t- e) = v(x, t) + e k v k (x, t) (7.9) 

fc>i 

with v(x,t) given by (|7.3J) satisfying (|7.1jl modulo 0(e 5 ). We will say that such a 
solution is monotone at the point x = Xq, t = to if 

u x (x , t ] 0) = v x (x , t ) ^ 0. 

According to the results of Section El all monotone solutions of the form (|7.9jl can be 
obtained by applying the transformation (|H.2jl to the nonperturbed solution ()7.2)1 (more 
precisely, one has to allow e-dependence of the function b(u)). 



Lemma 7.2 Let us perform the rescaling J7. 6] ) a/on^ wift 

e ^ A 7/6 e. 



(7.10) 



in i/ie quasitriviality transformation h3. ty) . Then the resulting solution to the perturbed 
PDE will be equal to 



u = v +\ l/3 + d 2 x 
(cf. H5.1}) ) where 



24 



c log v x + c e 



2 A I V xx 



^ XX ^ XXX ^ xxxx 

360 1920 vl + 1152^2 



c = c(v ), 

v = v(x,t) is the solution to the cubic equation J7.7| ). 



+0 (A 2 / 3 ) 
(7.11) 

(7.12) 
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Proof is straightforward. 

It remains to identify (j7.11j) with the formal asymptotic solution (|6.6J) to the ODE 
(jfi.lj) . This can be done by a direct substitution. An alternative way is to observe 
that, near the point of gradient catastrophe the perturbed PDE acquires an additional 
Galilean symmetry. Indeed, according to the previous lemma, locally one can replace 
the functions c(u), p{u) by constants c = c(v ), p = p(v ) (the constant p , however, 
does not enter in the leading term of the asymptotic expansion in powers of A 1 / 3 ). Let 
us show that in this situation any solution to the perturbed PDE of the form (|7.9j) 
satisfies also a fourth order ODE. 



Lemma 7.3 Let c(u) = cq, p(u) = po. Then for any solution u(x,t;e) of the form 
\7.9j) monotone at the point (xo,£o) there exists a formal series 

g(u; e) = g (u) + e k g k (u) 

k>l 

such that for arbitrary x, t sufficiently close to x , t the function u(x,t;e) satisfies, 
modulo 0(e 5 ), the following fourth order ODE 

x = t—l- + -f- 7.13 

ou{x) du{x) 

Here 

g'M = &(«)■ 



Proof It is easy to see that the flow 

u T = l-td x ^ (7.14) 

OU{X) 

is a symmetry of (|7.1jl . Combining this symmetry with one of the commuting flows 

u s + d x -^-r = 
du{x) 

one obtains another symmetry. The set of stationary points of this combination 

/ SHf 5H g , _ \ _ 
x \ Su(x) Su(x) J 

is therefore invariant for the t-flow. Considering the limit e — > it is easy to see that 
the integration constant vanishes on the solution (|3.2j) . (|7.2|) . The lemma is proved. 
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The ODE for the function u(x) is closely related to the so-called string equation 
known in matrix models and topological field theory (see, e.g., 0). Explicitly 

x = ta(u) + b(u) + c — {t [2a"u xx + a'"u 2 x ] + [2b"u xx + b'"u 2 x ] } (7.15) 



2 P0 (ta'> + b»)+±4(ta»' + b'") 



1 



4 P0 (ta"' + b"') + -^c 2 (ta IV + b IV ) 



U xxx Ua 



+ 



4p (ta IV + b IV ) 



l -p (ta v + b v ) + 



11 

1440 



1152 



cl{ta v + b v ) 
l(ta VI + b VI ) 



u xx u x 



u, 



Let us call the solution generic if, along with the condition k := —{a"'(vo)to + 
b'"(v )) 7^ it also satisfies 

c := c{v ) ^ 0. (7.16) 

Main Conjecture, Part 3. The generic solution described in the Main Conjec- 
ture, Part 1 can be extended up to t = t + 5 for sufficiently small positive 5 = 6(e); 
near the point (xq, to) it behaves in the following way 

, f e2c oY /7 rT f x ~ a o(t -t ) - X a' (t-t )\ ir) {4/7\ 

To arrive at the asymptotic formula 1)7.17)1 we do in (j7.15j) the rescaling of the form 
(|7.fi|) along with (|7.1()jl . After substitution to the equation (|7.15|l and division by A, 
one obtains 



x = a' ut 



y + 24 c ° + 2uu **) + j^'u" 



+ O (A 1 / 3 ) 



In derivation of this formula we use that the monomial of the form 

e U X U XX U XXX ' ' ■ 

after the rescaling will be multiplied by X D with 

D = \k + \{n + i 2 + ...) 
o 6 

due to the degree condition 

%x + 2i 2 + 3i 3 + ■ ■ ■ = k. 

Adding the terms of higher order k > 4 will not change the leading term. Choosing 

A = e 6 / 7 cf 

we arrive at the needed asymptotic formula. 

Clearly the above arguments require existence and uniqueness of the solution to 
(jl.8|) smooth on the real line described in the Main Conjecture, Part 2. 
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8 Concluding remarks 



We have presented arguments supporting the conjectural universality of critical be- 
haviour of solutions to generic Hamiltonian perturbations of a hyperbolic equation of 
the form (jl.2j) . In subsequent publications we will study the Main Conjecture in more 
details. The possibilities of using the idea of Universality in numerical algorithms to 
dealing with oscillatory behaviour of solutions to Hamiltonian PDEs will be explored. 
We will also proceed to the study of singularities of generic solutions to integrable 
Hamiltonian hyperbolic systems of conservation laws 



ut + d x 



dh{u) 



0. 



V 



det(r) ij ) ^ 0. 



U) 



Recall that, according to the results of [3T] the system (|8.1j) is integrable if it di- 
agonalizes in a system of curvilinear coordinates v k = v k (u), k = l,...,n for the 
Euclidean / pseudo- Euclidean metric 

n 

ds 2 = rjijdtfdul = g k (v)(dv k ) 2 , (77^) := (rj^y 1 , 

k=l 

v k t + \ k (v)v k = 0, k = l,...,n 

(in this formula no summation over repeated indices!). All Hamiltonian perturbations 
of the hyperbolic system (j8.1|) can be written in the form 



u\ + d x T] 



■• SH 

5uHx) 



0, H 



h{u) 



k>l 



e h h k (u; u x ,...,u 



dx, deghk 



k. 



We plan to study symmetries of the perturbed Hamiltonian hyperbolic systems. In 
particular, we will classify the perturbations preserving integrability and study the cor- 
respondence between the types of critical behaviour of the perturbed and unperturbed 
systems. The next step would be to extend our approach to Hamiltonian perturbations 
of spatially multidimensional hyperbolic systems (cf. [H]). 



Appendix: Bihamiltonian structures associated with 
the perturbations of the Riemann wave hierarchy 

Theorem A.l For arbitrary two functions c = c(u), q = q(u) the family of Hamilto- 
nians \2. .5)) with 



p(u) 



960 



c q' 



s(u) = 0. 



(A.l) 



is commutative 



{H f ,H g } 1>2 = (mod 0(e 6 )) V/ = f(u), = g( 



u 



(A.2) 
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with respect to the Poisson pencil of the form 

{u(x),u(y)} 1 = 5'(x - y), 

{u(x) , u{y) } 2 = {u(x) , u{y) } ™ + e 2 {u{x) , u{y ) } M + e 4 { u (x) , u(y) } [4] + 0(e 6 ). 
i/ere i/ie terms of order 0: 

{u(x), u(y)} [ 2 ] = q(u)S'(x - y) + -q'(u)uj(x - y) 

All terms of higher orders are uniquely determined from the bicommutativity 
provided validity of the constraint hA.l\) . Namely, the terms of order 2: 

{u(x), u(y)}f = ^-S'"(x - y) + ^- (cq')'uj"(x - y) 
8 16 



+ 



+ 



c"q' + c'q" + 5cq" 



16 



6 



48 



2 i c q u xx 



c"q" + c'q'" 



cq 



(4) 



48 24 
The terms of order J^: 



48 



u x + 
1 



16 



+ 



7cq"u x 
48 



u x +-(cq +cq )u x u xx + —u xxx 



S'(x - y) 

5(x - y) 



{u(x), u(y)}f = {SccW + c 2 q") S V (x-y) + ^ {3cc'q' + c 2 q")' u x 5 IV (x - y) 



+ 



+ 



Vc'V cc^V 3c' 2 q" 5cc"q" cc'q" 2 u x 2 c 2 q" 3 Wcc'q'" 3c 2 q"q"' , c 2 q^' 



32 



32 



32 



IS 



240g' 



+ 



VV 3ccV , 17«V cY' 2 



64 



64 



+ 



192 



480g' 



+ 



480g /2 

19c 2 g w 
960 



192 



u. 



640g' 
5"'{x-y) 



64 



W c'c"'q' cc^q' Wc'd'q" 23cc"'q" hcc'q^ 7cc"q"' c 2 q^ 3c' 2 q"' 



128 



32 



128 



128 



384 



64 



64 



96 



32 



c' 2 q" 2 cc"q" 2 cc'q" 3 c 2 q" A I7cc'q"q"' 21c 2 q" 2 q"' 9c 2 q"' 2 9c 2 q"q^' 

+ — .o - ~ ,3 7T77. ; H 



160g' 160g' 80g' 2 160g 



640g' 



1280g 



/2 



1280g' 1280g' 



+ 



9c'c"q' 3cc"'q' llc'Y' 13cc"q" 3cd 
+ „, + — ~~~ + 



64 



6 I 



64 



64 



3c 2 q' 

+ ... .o + 



69ccY" 13cW" 3c 2 a (4) 



160g' 320g 



/2 



320 



640g' 



+ 



80 



U T . Unt 



/ c' q' cc"q' I3cc'q" c 2 q"~ era 



32 



32 



192 



320g' 



60 



o"(*-y) 



+ 



<•""(/" r'c'"(i" (■(-■'"( I" c'd'q" 2 



48 



32 



96 



cc"'q" 2 c' 2 q" 3 cc"q" 3 



160g' 480g' 160g' 2 160g 



/2 



/ //4 
CCO CO 

' + 



80g 



/3 



160g 



+ 



35c'c"g'' 



384 



hod" J" 9c' 2 q"q"' 9cc"q"q"' llcc'q" 2 q"' 13c 2 q" 3 q"' cc'q'" 2 19c 2 q"q"' 2 Yld 2 q^ 



128 



640g' 



640g' 



320g 



/2 



640g 



/3 



64g' 1280g 



/2 



384 



20 



5cc"q^ cc'q"q^ 17c 2 q" 2 q^ UcVVg 35cc'q^ 11c 2 q"q& c 2 q^ 



96 



64g' 



+ 



128 



+ 



32 



+ 



128 



1920g' 2 
+ 



1280g' 



1152 



3840g' 288 



3c" V c'c"'q' cc( 4 )(7' 91cVY' 37cc , y c' 2 q" 2 



384 



+ 



384 



// 1/2 / //3 

cc f cc q 



60q' 60q' 30q' 2 60q 



c 2 q" 4 59c' 2 q"' 



/3 



320 



53cc'Y" A7cc'q"q"' 173c 2 q" 2 q"' 77c 2 q"' 2 169ccY 4 ) 77c 2 q"q^ 73c 2 q^ 




+ 



3840g 



/2 



cc"Y 5c'Y' cc'Y' 



3840g' 



/ //2 

cc g 



+ 



cY' 3 



+ 



960 3840g' 2880 

157cc'q"' 5c 2 q"q"' 31cV 4) 



Ux U xx 



+ 



+ 



96 ' 16 80g' 160g' 2 ' 1920 384g' ' 1920 

l2 q" 3cc"q" cc'q" 2 c 2 q" 3 19cc'q'" llc 2 q"q"' llcY 4 )\ 



12 



+ 



, c' 2 q> cc"q' 
\ 128 128 



32 60g' 
llcc'q" 



+ 



120q 



,2 



+ 



160 



c 2 q" 2 



384 



320g' 



+ 



17c 2 q'" 
1920 



640g' 
S'(x-y) 



+ 



480 



) 



tlx Ua 



c" 2 q"' 




c'c"'q"' ccWq'" 
+ 384 

//„/// 2 



CC 



+ 



_ c'c"q"q'" 
~ 640g' 

3cc'q"q"' 2 



cc"'q"q"' c' 2 q" 2 q"' cc"q" 2 q"' 



1920g' 



640g 



/2 



640<? 



i2 



ccVV 

/3 



320g 



640g' 640g' ' 640g' 2 

ccq"~q^' c"q"~q( 



c 2 q" 2 q"' 2 c 2 q"' 3 7c'c"q^ cc"'q^ 

+ .... .„ + _ + 



,/2 » ( 4 ) 



320g 



/3 



1280g 



/2 



„//„//„(4) r/ V 2 „(4) r .2„//3 n (4) 3cc 'g'" g (4) 



640g' 



320g 



/2 



640g* 



/3 



640g' 



ljcVYVj c 2 g w 



3840g' 



/2 



384 128 640g' 

17c' 2 q^ 
1280g' ' 2304 



5ccV 5 ) cc'q"q^ cY'Y 5) c 2 q"'q^ 5ccY 6) c 2 q"q^ c 2 q^' 



576 



640g' 1280g 



,/2 



960g' 



1152 



3840g' 2304 



+ 



V'Y cV'V' cc( 4 )g" c'c'Y' 2 cc w g" 2 c'Y' 3 cc'V' 3 

H ~ I ; 7777. , I" ,2 + 



CC 



/„'/4 



64 



48 



192 



160?' 480g' 160g' 



160tf 



/2 



80g 



/3 



+ 



160g 



/4 



97cVV; 13cc"Y" AT cc"q"q"' 19cc'q" 2 q"' 11c 2 q"\"' cc'q'" 2 3c 2 q"q"' 2 

+ ' + 7777^, 7^-=- + —7 ^-^+ ' ' 



960 



320 



60g' 



60g' 



480<?' 



/2 



480g* 



f 3 



A8q> 



160q 



.12 



19c'V 4 ) 67cc"^( 4 ) cc'q"qW llc 2 q" 2 q^ c 2 q"'q^ 131cc'q^ c 2 q"q^ c 2 q^ 




A8q' 



n /2 1,2 

7c q 



+ 



960g 



7 2 



80g' 



+ 



2880 



+ 



384 960<?' 
13c 2 q'"q"' 3c z q 



960g' 480g /2 960g /3 



7c 2 q" A 59c' 2 q"' 



960 



240g' 

23cc"g w 
320 



+ 



180 



cc'q"q"' 
30q' 



2„„2jn Q „ 2 „///2 131cc / ? (4) 3c 2 q » q (4) 3 l c 2 g (5) 



640g 



,2 



320q' 



+ 



1920 



320g' 



+ 



2880 



+ 



3c'c"q" cc"'q" 



64 



64 



,2 n2 it i/2 i //3 

c q cc q _^ccq 



160?' 160?' 80g 



i2 



c 2 q" 4 47c' 2 q"' 13cc"q" 



160g- 



/3 



960 



240 



21 



13cc'q"q'" c 2 q" 2 q"' 7c 2 q"' 2 A9cc'q^ 7c 2 q"q^ 23cVM 2 
A80q' + 60q' 2 960g' + 960 960q' + 2880 ) % 

/ 5gV 5ccV_ccV^ cV' 3 3ccV^ _ cVV^ c 2 g (4) \ 
+ 1 192 + 192 96g' + 192g' 2 + 64 96</ + 96 ) UxxUxxx 

( c' 2 g" cc"q" cdq" 2 c 2 q" 3 9cc'q'" c 2 q"q w c 2 q^ \ 
+ I ~64~ + 64 ~ T607 + 320g' 2 + ~320 160g' + 160 ) UxUxxxx 



+ 1 192 960g' + 480 J 



5(x - y) 



To prove the Theorem one has to analyze the commutativity conditions 



Su(x) 5u(x 

for arbitrary two functions f(u), g{u). Here 

1 e 2 
L = qd x + -q'u x - —cq'dl + ... 

is the Hamiltonian differential operator associated with the second Hamiltonian struc- 
ture. To prove validity of Jacobi identity one has to check that the e-terms in the 
second Hamiltonian structure can be eliminated by the quasitriviality transformation 
described in Section |3J We will omit the calculations. 

Observe that the family of bihamiltonian structures given in Theorem I A . 1 1 dep ends 
on two arbitrary functions c = c(u), q = q(u), in agreement with the results of j2H|. It 
is understood that the Jacobi identity for the Poisson pencil holds true identically in 
A modulo terms of the order 0(e 6 ). 
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